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We investigate the free energy and the penetration depth of a quasi-two-dimensional d-wave 
superconductor in the presence of a weak magnetic field by taking account of thermal, nonlocal, 
and nonlinear effects. In an approximation in which the superfluid velocity v s is assumed to be 
slowly varying, the free energy is calculated and compared with available results in several limiting 
. cases. It is shown that either nonlocal or nonlinear effects may cut off the linear- T dependence of 

' both the free energy and the penetration depth in all the experimental geometries. At extremely 

low T, the nonlocal effects will also generically modify the linear H dependence of the penetration 
depth ( "nonlinear Meissner effect" ) in most experimental geometries, but for supercurrents oriented 
along the nodal directions, the effect may be recovered. We compare our predictions with existing 
^ , experiments on the cuprate superconductors. 

^ PACS Numbers: 74.25.Nf, 74.20.Fg 

s : 

^ ■ I. INTRODUCTION 

^— > . 

Ctf . 

Many unusual aspects of the unconventional superconductivity in the cuprate and heavy fermion materials are 
related to the existence of nodes in the energy gap. Although intensely discussed in the 1980s and 1990s, most of 
these phenomena were already known much earlier to students of the unconventional superfluid 3 He, which exhibits 
several phases, some with gap nodes and accompanying low-energy fermionic excitations. For example, power laws 
in specific heat, transport properties, and NMR were predicted and measured in the anisotropic, high-pressure "A" 
phase, along with singular responses to impurities, unusual vortex phases, and even anomalous Josephson currents 
jij. From this perspective, the one important difference between the new superconductors and the 3 He system is that 
they are charged; unconventional superconductivity can be uniquely probed by studying the response of the system 
' to an applied electromagnetic field. 

Even the simplest such phenomenon, the expulsion of a weak applied magnetic field, or Meissner effect, is still the 
subject of intense discussion today in the context of unconventional superconductivity. The new materials are strong 
type-II superconductors, with large values of the London penetration depth A and small coherence lengths t;o, such 
that the condition Ao ^> £o is expected to hold. It was therefore anticipated in the earliest theories of the Meissner 
Q\ \ effect in unconventional superconductors Q that the electrodynamics of the system could be treated as local, i.e., 
as though the Cooper pairs were point objects, neglecting the spatial variation of the electromagnetic wave over the 
extent of the pair. Under this rather plausible assumption, the simple anisotropic extension of the BCS theory is 
expected to hold, 



o 
o 



o 



where the quasiparticle spectrum is Ek — ->J£f. + Aj£, is the normal-state electron band, is the momentum- 
dependent order parameter, / is the Fermi function, k\\ is the direction of the supercurrent, and (• • -)fs represents a 
Fermi surface average. The London penetration depth at T = is given by Ao = \J mc 2 / {Aitne 2 ) , and is typically one 
to several thousand Angstrom in the new materials. 

At low temperatures, the local penetration depth calculated from Eq. (1) may be shown Q to vary as a power law 
of temperature, A(T) — Ao ~ (T/Ao) Q , where Ao is the gap maximum and the exponent a depends on dimensionality, 
nodal topology, and the rate at which the gap goes to zero near a nodal wave vector k n . For line nodes in a three- 
dimensional (3D) system (or point nodes in a 2D system) like the d x 2_ y 2 state thought to characterize the cuprate 
superconductors, a is found to be 1. Thus when Hardy et al. H measured a linear temperature dependence in the 
YBa2Cu306.95 (YBCO) system down to a few degrees K, it was dramatic evidence in favor of line order parameter 
nodes and possible unconventional pairing. 
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Recently, Kosztin and Leggett M questioned the theoretical basis for this identification, pointing out that the 
approximation of local electrodynamics used to derive Eq. (1) is not valid for the quasiparticle states near the 
nodes whose occupation determines the asymptotic low-temperature behavior. The coherence length £o = vf/ttAq 
determines the scale of order-parameter variations, as in the usual BCS theory, but an effective coherence length 
£ofc = vp/irAk appears in the electromagnetic response, as noted in Ref. |5|. The divergence of £ofc at the nodes means 
that quasiparticle states near the nodes must always be treated nonlocally ||. A full nonlocal calculation gives a 
crossover from a SX ~ T regime for E non i oc <C T -C Ao to a SX ~ T 2 for T <C E non i oc . In the YBCO system, this 
crossover scale E non i oc ~ Aq^q/Xq is about 1 K, and thus the identification of the d-wave state on the basis of the 
local theory is not really in question; measurements on the purest YBCO crystals down to 1 K or so have seen no 
deviation from the linear behavior. JtJ Quadratic dependences 6X ~ T 2 measured in some samples are more likely due 
to impurity effects |^||[|, which also produce aT->T 2 crossover. 

An interesting perspective on the question of the linear T dependence of the penetration depth was provided 
recently by Schopohl and Dolgov JIc| l , who observed that if the d-wave penetration depth temperature dependence 
were to remain linear down to T = 0, one would be unable to satisfy the third law of thermodynamics (Nernst's 
theorem). In the framework of linear-response theory, they expressed the change of the free energy in the presence of 
the magnetic field (for fixed external current source) as a functional of the penetration depth, 



AT (SD) = _J_ f A 
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where A q (T) is the Fourier component of the magnetic vector potential. The resulting entropy is 

d 3 q d 



) 3 dT 



A 2 (q,T) 



|A q (T^. (3) 



The Nernst theorem requires S to be when T — > 0, inconsistent with A(T) ~ T. 

Schopohl and Dolgov speculated that this result might imply the instability of a pure d-wave state within BCS 
theory at sufficiently low temperatures. It was pointed out in Ref. [ JLj] , however, that the extra magnetic-field-induced 
quasiparticles that deplete the shielding current, leading to a magnetic- field-dependent penetration depth known as the 
nonlinear Meissner effect |l2],[l3| , were neglected in this analysis. Several authors have recently considered the effects 
of these nonlinear terms in the supercurrent on transverse magnetization measurements fl4|| and on the structure of 



the vortex lattice p5 16|. The basic idea is that, like impurity effects, the magnetic field itself may serve as a Cooper 
pair breaker that creates nodal quasiparticles, leading to a temperature dependent penetration depth SX(T) ~ T 2 at 
temperatures below the scale for nonlinear electrodynamics, E non u n ~ v s Icf with v s a typical supercurrent velocity and 
hp the Fermi wave vector. The current authors p7| then provided a synthesis of the various arguments given above, 
pointing out that for any experimental geometry with given disorder, external current distribution, and Ginzburg- 
Landau parameter n, a sample of d-wave superconductor would inevitably avoid the violation of the Nernst theorem 
by creating a T 2 term in its penetration depth through a competition of nonlinear, nonlocal, and impurity effects. 

Although this issue of principle has been resolved, an important prediction of the d-wave model for the electro- 
dynamics of high-T c superconductors (HTS) has not yet been confirmed. In the early discussion of the symmetry of 
the HTS order parameter, Yip and Sauls proposed that the angular position of the gap nodes could be probed by 
a measurement of the magnetic field dependence of the penetration depth. In the local limit and for T — > 0, they 
predicted that the nonlinear effect induces a linear-H term in the penetration depth, 

SX(H) _ X(H)-X(0) _ f H/H v s || node 

A(0) ~ A(0) ~ { ^ H / H o v s || antinode ' I > 

where H is the applied magnetic field at the surface and 

is of order the thermodynamic critical field of the system with <!>o the flux quantum. The prefactor of the linear 
H term in Eq. (^) depends generally on the angle that the supercurrent makes with the crystalline d axis and is 
independent of T (the T dependence appears in the H 2 term as discussed recently by Dahm and Scalapino [Q ) . 
Several experimental groups tried to verify the prediction of Yip and Sauls, but failed to identify a linear H term 
which scaled with the temperature according to the theory |I8|-pO| . 
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There are three relevant energy scales in the low-energy sector in the Meissner state, T, E non i in , and E non \ oc . The 
first two are experimentally controlled parameters while the last is an intrinsic one. Obviously, A(T) ~ T behavior 
should be expected at T ^ E non i oc ,E non u n . But in the opposite limit, i.e., at extremely low T, either nonlinear or 
nonlocal effects may play a crucial role in modifying the linear T behavior. Thus a full study including both nonlinear 
and nonlocal effects in the Meissner state is necessary. A similar investigation has been performed by Amin, Affleck, 
and Franz |l5f| in the mixed state where H > H c \, the lower critical field. It was found in this numerical study that 
the nonlinear and nonlocal corrections were of the same order. This result is not surprising: in the mixed state the 
nonlinear energy scale and the nonlocal one are not independent, since typical spatial variations take place on a scale 
of the magnetic length or intervortex distance, fixed by the fluxoid quantization. 

In this paper, we perform a full study of a quasi-two-dimensional (2D) d-wave superconductor in the presence of 
a weak magnetic field. We start from a phcnomcnological BCS Hamiltonian with a quasiparticle energy spectrum in 
the normal state, £k = et — /i (p is the chemical potential) and a d x 2_ y 2-wave order parameter, A v — Aq cos2i^ (ip is 
the angle of k with the a axis). Using the functional- integral approach, we calculate the free energy and compare it 
with that obtained in other theories. Our main approximation is that the supercurrent along the surface boundary is 
a slowly varying function of the distance to the surface. Thus we take v s = (e/mc)A q= o as the supercurrent to which 
the quasiparticles couple, leading to the Doppler shift in the energy spectrum of particles and holes ±E^ + k ■ v s , 
where E k = + A*] 1 / 2 . The response of the quasiparticles to the A q ^o modes is treated perturbatively, and the 
resulting response function involves both the nonlinear and nonlocal effects. From this response function we then 
calculate the penetration depth, and express it as a scaling function of the two parameters E non i oc /T, E non y m /T. We 
show explicitly that, below a temperature T*, the linear T behavior of the penetration depth is always modified 
by either nonlinear or nonlocal effects to a quadratic function of T, while the linear H behavior in the penetration 
depth is also cut off by the nonlocal effects when H is below some characteristic critical field H* , which is geometry 
dependent. For most geometries H* ~ H c i, implying the unobservability of the nonlinear Meissner effect, as discussed 
in our recent paper on this subject p2| . But for some special geometries, in which the supercurrent is along nodal 
directions, the nonlinear Meissner effect may be recovered. We also present a brief discussion on the geometry in 
which the magnetic field is along the conduction plane, which presents some special features due to the highly layered 
nature of the cuprates. Finally, we investigate the consequence of the spatially varying nature of the supercurrent 
and find that this property does not influence the qualitative behavior of the penetration depth. 

The paper is organized as follows. In Sec. II, the model Hamiltonian and the partition function are obtained. In 
Sec. Ill, the expression of the free energy as well as its comparison with various theories in several limiting cases is 
presented, and the supercurrent and the response function are discussed. Detailed study of the penetration depth in 
the presence of a constant v s is given in Sec. IV. Section V is devoted to concluding remarks and comparison with 
existing experiments. Finally, a discussion of the renormalization of the penetration depth due to space-dependent v s 
is included in the Appendix. 



II. MODEL HAMILTONIAN AND PARTITION FUNCTION 



We start from a phenomenological d x 2_ J/ 2-wave BCS mean-field Hamiltonian of a 2D system in the presence of a 
weak magnetic field p3j: 
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iV--A(r) -ii}c a (r) + I d 2 rd 2 v' A(r, r')c{(r)c{ (r') + H.c 



-|A D 



+ J d 2 rd 2 r'V(r-r')\b(r,r' 



(6) 



where — V(r — r') < is the effective attraction between electrons responsible for the d x 2_ y 2 pairing, A(r, r') = 
— V(r — r')fe(r,r') is the d x 2 _ y 2 -wave pairing order parameter with 6(r,r') = (c^(r')c|(r)), and A(r) the magnetic 
vector potential. It is well known that in the Meissner state the supercurrent formed by the Cooper pairs effectively 
screens the magnetic field from penetrating into the bulk of the sample. Within the penetration depth away from 
the surface, both magnetic field and the supercurrent decay with the distance to the surface, and the pairing order 
parameter acquires, in principle, a spatially varying phase. As usual we write A(r,r') as A(r — r')e l ^ B ^ with a real 
A(r — r') and R = (r + r')/2. Note we neglect the supression of the magnitude A at surface of the sample, since we 
expect these corrections to the penetration depth will be small in the limit Ao » £0 of interest. 
The partition function of Hamiltonian (pf) is 
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where 



Zq = exp 



N S A 2 

Vn 



(7) 



(8) 



with N s the total number of superconducting electrons, rfi(r) = Cj(r), ^(r) = c|(r), and VW a /3(r, r') is a matrix 
element of the 2x2 matrix M in the Nambu representation, 



M(r,r') 
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i^ + §A(r) -h\5[y- r'i 
A(r - r ') e -^( R ) 

The use of a unitary transformation 

C/(r,r') = 



A(r-r')e 



(9) 
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e #(R)/2 q 

e -#(R)/2 



(10) 



to A-i is made for eliminating the phase factor of A(r,r'), leading to 
M(r, r') = U'< (r, r')M(r, r')U(r, r') 



V 

where 
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A(r-r') 
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iV ■ v s (r) - |mw?(r)| 5(r' - r) 



, (11) 



v.(r) 



eA(r) V0(r) 
mc 2m 



(12) 



is the supercurrent velocity. In writing down Eq. (fil"l ) we have used the relation V • v s (r) =0 which corresponds to 
the conservation of the supercurrent. v s (r) should be gauge invariant, which is guaranteed by the compensation of 
the two terms in the right-hand side of Eq. (|l2| ) under a gauge transformation. It is necessary to fix the gauge before 
proceeding. In the Meissner state the most convenient gauge choice is the London gauge in which V0(r) = 0, and 
hence 



v s (r) 



eA(r) 



mc 



V • A(r) = 0. 



The partition function Z in Eq. (0) can be expressed in momentum space as 



Z = Zq exp 



£tr In (Q^(iuj n ) + Q 



(2: 



where uj n is the fermion Matsubara frequency. 



-iw„ + Ck + M k}k A k 

A k -uo n - Ck - M_ k ,-k 
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-M_ k '-k 



(13) 
(14) 

(15) 
(16) 



with Ak = J d 2 rA(r) exp(ik-r) ~ A cos 2ip on the Fermi surface, M kl ^ 2 = — (e/mc)k 2 • A kl -k 2 + (e 2 /2mc 2 ) A p • 
Ak!-k 2 -p with Ak the Fourier component of 5A(r) = A(r) — A (A is the spatial average), £k = £k — A* i s the energy 
spectrum of an electron in the normal state, and Q{iuj n , k) is Green's function matrix in Nambu's representation 



4 



(/(wj„,k) 



At 



A k 

-£k--fl/k.k 
Wk 



Here 



where 



W k = -(iw n + k-v s ) 2 + ^, 



eA 



(17) 



(18) 



(19) 



and £ k = [(£k + (e 2 /2mc 2 ) / (d 2 p/(27r) 2 )A 2 ) 2 + A 2 .] 1 / 2 . In Eq. (|||), we have explicitly separated the contributions 
to the free e nerg y from the homogeneous (q=0) superflow and those corresponding to q ^ components. 
Z in Eq. dlj) can be written as 



Z = Z cxp < tr lnQ*- 1 -* (iui n ) > exp < tr In 1 — G(iuj n )Q 



)(2) 



2q ■ Z\ ■ Z 2 , 



(20) 
(21) 



Zx = exp{^tr lnQW(iw n )} = J]{1 + e -(k-v.+^)/r } . + e -(k.v„-^)/T } _ 

n k k 

Z 2 may now be calculated by expanding in powers of A q ^o • To second order, we have 

hL§2 = J2 tT - G(^n)Q (2) ) ~J2 tT -S(^n)Q (2) - \ (-G{lOJ n )Q {2) ) 2 

= -^L J (^j2 y (2^)2 {en(iw„,k + )0 u (^„,k_)(k_-A q )(k + -A_ q ) 

+&2(io; n , k + )022(iw n ,k_)(k + ■ A q )(k_ • A_ q ) + Q\i{iw n , k+)£ 2 i(«w n , k_)(k+ ■ A q )(k+ • A_ q ) 
+g 12 (iw n ,k + )g 2 i(w„,k_)(k_ ■ A q )(k_ ■ A_ q )}, (22) 

where k± = k ± q/2. Note that only the paired electrons, i.e., the electrons near the Fermi surface, can produce a 
diamagnetic effect. Thus we can replace k that couples to A q and v s by kp. Adopting the London gauge q • A q = 0, 
we simplify lnZ 2 in Eq. rt22) as 



lnZ 2 



E 



d 2 k f d 2 q {iuj n + k F ■ v s ) 2 + + A+A_ 



(2tt) 2 J (2tt) 2 
where £± = £k± , A± = A k± and W± = W k± . 



(k F • A q )(k F • A_ q ), 



(23) 



III. FREE ENERGY, SUPERCURRENT AND RESPONSE FUNCTION 

The total free energy density of the superconductor in the Meissner state T is now given by the sum of contributions 
from both electrons and magnetic field, where the former can be calculated in terms of the partition function. We 
obtain 



where 



T=-T\nZ + ^- / d 2 r|Vx A(r)| 2 ^^+J- 2 , 

07T 



T x = -nmv 2 s -T ( (in [l + e -M.+W] + In [l + e -i^-E k )/T\ j 
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with 
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rfn v T) — ^ (a 4- 2T V [ ^ (V ^n+^F ■ V s ) 2 A + A _ \ 

/C(q ' Vs ' T) -^( v 1 + ^^ J / (2^F (kF) H W+T^ J • (27) 

^-"i describes the kinetic energy of the quasiparticles with energy spectrum ±-E& shifted by k F • v s . Formally, it is 
nothing but the Doppler shift of the quasiparticle with proper momentum k in the lab frame with local superfluid 
velocity v s (r). @ 

It is interesting to look at T2 in Eq. (|2^). In the linear limit v s — > 0, /C(q, v s ,T) in Eq. (|27j ) reduces to the usual 
linear response function of a superconductor in the Meissner state p5[ . The fact that /C(q, v s , T) is indeed a response 
function can be seen from minimization of the free energy with respect to A q , dTijdK^ = 0, which gives the Fourier 
component of the current inside the superconductor, 

jsc(q) = -£(q,v a ,T)A q , for q ^ 0. (28) 

Therefore Eqs. (p4|)-(p7|) describe the Bogoliubov quasiparticles with the Doppler-shifted energy spectrum responding 
to a weakly spatially varying magnetic field. Note that /C(q, v s ,T) is obtained in an infinite sample. In the presence 
of a surface boundary, an external current j cx t(q) is introduced to coincide with the real boundary condition (see 
discussion in Sec. IV), leading to the total current as 

jtot(q) = ^9 2 A q - j ex t(q) - £(q, v S) T)A q . (29) 



a-b plane 




FIG. 1. Geometry of (010) surface case. 
By summing over the Matsubara frequencies, Eq. (|27j) becomes 
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/C(q,v s ,T) = — ^ + J/C(q,v s ,T), (30) 
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[/(^-k F -v s ) + /(^ + k F .v s )]A 2 
^ 2 - A£{A 2 + [(q • k F )/(2m)] 2 - 

where f(x) is the Fermi function. <5/C(q, v s ,T) is geometry dependent. The first two terms in Eq. (32) represent 



+2 {k F )\ / d^Re 7 k (32) 

' V^A 2 :{A 2 +[(q.k F )/(2m)P-^}/ FS / 
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the nonlocal renormalization of the T = response, while the third represents the combined nonlocal and nonlinear 
corrections to the T dependence. For a real quasi-2D system, like the high-T c compounds, with the ab plane as the 
conduction plane, we first consider the geometry where the magnetic field is parallel to the c axis and thus v s and 
the penetration direction q are in the ab plane, and in general, v s makes an angle 9 with the a axis. This geometry 



G 



is shown in Fig. 1. Note that for v s /\\ antinode the possibility of the formation of an Andreev bound state exists. 
Walter et al. (2^| have argued that such states can have significant effects on the low-T penetration depth for v s | 
node. The effects on the field dependence of the penetration depth are beyond the scope of the current formalism, so 
results for these directions should be treated with caution. We will consider a different situation, where the magnetic 
field is in the ab plane, at the end of Sec. IV. 

An examination of Eq. ( |32|) shows that the small regions around the four nodes give the main contribution when 
v s kF,T <C Ao. This reflects the low-energy quasiparticle excitations near nodes. It is straightforward to integrate 
out the angle of the Fermi wave vector, leading to 



6tC(q,v s ,T) 



c T ^ 2 au e t eugi 



fcrA§A ^ 



T 



T 



(33) 



where a = qJtF /(2-\/2m), e — v s kp/V2, u$i = | cosf? + Zsin#|, and F\(z\, z^) is a two-parameter scaling function 

1/2 



F\(zi, z 2 ) = - Zl + [ln(e 22 + 1) + ln^ 2 + 1)] 



dx 



1 



I 



1 e x+Z2 + 1 



1-1^ 



(34) 



Now we may study two limiting cases. 



A. v 3 k F < T 



In this case we can expand T\ in Eq. (25) up to the leading order of v s kp/T 



<ri — r BCS 



1 2 If ePk 



2 nmV ° T J (Jhr)» 



(I + e <VT )2 _ 



(35) 



where E { ° ] = (£f + Af) 1 / 2 and ^ s is the free energy density of the superconductor in the absence of a magnetic 
field, 



tt(°) - _t 

-^BCS — 1 



In(l + e- E f ) ' T ) + ln(l + e E ^ l T ) 



Inserting Eq. ( |l9| ) into Eq. fl35| ) yields 



•^i^^Bcs + ^(0,0,r)A 2 , 



with £(0,0, T) 



4ttA,t 



5tC(0,(i,T) and 



nA| 
V 



(36) 



(37) 



(38) 



As for T2 in Eq. (|2£]), up to the leading order of A q , the response function /C(q, v s , T) is simply replaced by /C(q, 0, T), 
the usual nonlocal response function It is straightforward to see from Eqs. ( |33| ) and ( p4"| ) that 



<5/C(q,0,T) = <5/C(0,0,T) 



8ttA^ 



E 



u ei \ ri"' ' - IT f dxf(x) 
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aue-i €T, J = ±1 
aug 1 > T, £ = ±1 



(39) 
(40) 



a special case of which, 6* = 0, has been presented in Ref. M. In Eq. ( |4C| ) £(3) ~ 1.20 is Riemann's zeta function. 
Therefore the change of the free energy due to the presence of the magnetic field is 



8tt J (2tt 



(41) 
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Comparing Eq. ( j4l| ) with Eq. (g), we seem to find a sign discrepancy. However, notice that the free energy in Eq. (g) 
is calculated at a fixed external current and hence is the Gibbs free energy, whereas that in Eq. ( f4l|) is the Helmholtz 
free energy, A.F( Holm ), since a fixed A q is assumed in the calculation. It can be checked that these two kinds of free 
energy are related to each other through the Legendre transformation 

A ^( SD ) = A ^(H clm) _ I J |^j ext (q)-A q = -i- J ^b 2 + ^)( 9 ,0,T)]A q . (42) 



In the last equality, Eq. (|29|) was used. 

The entropy can be calculated from the Helmholtz free energy in Eq. ( pi] ) with A fixed, which is easily shown to 
be identical with that in Eq. (m 



S(T) = - 



dT 



= S {SD '>(T). (43) 

fixed A 



Equations ([II]) and ( |38"1 ) imply that in the local limit a — > 0, AJ 7 oc T. If this linear T behavior were to hold at 
T — > 0, the third law of thermodynamics, the Nernst's theorem, would be violated as pointed out by Schopohl and 
Dolgov pj. However, the intrinsic nonlocal effects may cut off the linear T term [(||l7|]. Kosztin and Leggett showed, 
in the special geometry where 9 = 0, that the linear T behavior of SIC(q, v s — > 0, T) will be changed to a higher power 
of T when T is smaller than a. In the case of more general 9, the linear T term is preserved until T is smaller than 
both auei,l = ±1. For v s along any node, u$i = for I = 1 or — 1. The nonlocal effects thus disappear along this 
nodal direction. In this case, the linear T singularity may be cut off by another effect, the nonlinear effect, which we 
now discuss. 



B. Vsk F » T 

Notice that in Eq. ((3^), s is always accompanied by uei, which, for 9 near any node, may be very small. We 
consider first the situation where 9 is not close to any node, so that eugi S> T. The small perturbation parameter is 
now T/(v s kFUgi). The kinetic free energy T\ becomes 

1 f°° f <i 2 k 

T x = -nmv 2 s - 2T / dxN(x)ln[l + e -^ T ] - / -—j E k , (44) 

2 J-oo J ( 27r ) 

where 

N[x) = f j!r^5[x ~ (v. • k F + E k )} = -^r-J2( x + £u eiW x + £u ez), (45) 

with i9(y) the step function and N the density of states on the Fermi surface in the normal state. Inserting Eq. ( |45| ) 
into Eq. (f|J) we get 

T 77(0) , 1 2 1 Ce nmk F 3 tt 2 nm \^ 2 

Tx = ^ BCS + tT ' ~ 3^W" S " IAo^f ^ ' (46) 

where 

Ce = - ^ i4 = - (|cos(9 + sm(9| 3 + | cos 9 - sin (9| 3 ) . (47) 
l=±i 

The qualitative behavior of T\ as a function of both, v s and T, is consistent with that obtained by Volovik [pr] . 
Here we see that the prefactors of both the v$ and v s T 2 terms are 9 dependent. 

The scaling function F\{z\, z%) in Eq. (|i]) has the following asymptotic behavior at Z2 ^> 1: 

( Z2 Z\ < Z 2 

F A (zi,Z 2 )^ „ S1 i » 3 («3-Hr a ) . „ ? • (48) 
T + 6^ ^2 < Zl • 



This may be inserted into ( p3[ ) to find that <5/C(q, v s ,T) = X)z=±i ^^Hli v s j ^) > where 
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It is clear from Eqs. ( f40| ) and (|49]), that the linear T term in 5/C(q, v s , T), and hence in J^i is modified to a higher power 
of T by either nonlinear or nonlocal effects. Combining this result with T\ in Eq. (44), we see that a thermodynamic 
instability, associated with the violation of the Nernst theorem, is avoided. 

For 6 near a nodal value, ugi ~ for I = 1 or — 1. This means that those quasiparticles in nodal regions with 
wavevectors nearly perpendicular to v s contribute negligibly to the nonlinear effect such that the linear T singularity 
in the free energy may not be cut off. These quasiparticles, however, acquire large nonlocal effects since ql v s , 
which may effectively modify the T dependence of the free energy. Similarly, it follows that the quasiparticles which 
generate a small nonlocal effect produce much larger nonlinear effects. 

With the expression for the free energy, we can easily obtain the supercurrent. The part of the supercurrent arising 
from the kinetic free energy T\ in Eq. ([46[) is given by 



• (kin) 



e dT\ 



1 d( e 



2A Vv^ 3\/2 88 



where z || c. It is clear that due to the nonlinear correction jic*™' 1 is not parallel to v s except for v s || node or antinode 



z x v. 



(50) 



in the case of which we recover the Yip-Sauls nonlinear supercurrent form ]l^ , |T3| ], 



• (kin) 



-env a (l - w „ 



v.. 



node 
antinode 



(51) 



The fact that jsc mj is not parallel to v s for general 9 leads to the interesting magnetic torque phenomenon that has 
been discussed in Ref. Q. 

j(q) can be obtained from Eqs. (p8|), (E 



j(q->0) = 



and ( f49|) . It is nontrivial to examine the supercurrent in the local limit 

(52) 



»0j 



■4%(1-^) A < 
c n _ 1 v„k F \ a 



v.. 



node 
antinode 



If we simply take (e/mc)A q ^ as v s , we find that j(q — > 0) is of the same structure as in jic*"' in Eq. (|5l|), but with 
an extra prefactor 2 in the nonlinear correction term. Note that j(q = 0) in the present theory is not equivalent to 

jsc > but instead, 



•1 m e dT 

j q = = — t— 

m ov s 



'(kin) 



2mc J (2tt) 2 

If /C(q, v s ,T) can be replaced by its local-limit value, Eq. (|53|) becomes 



d 2 q 5/C(q,v s ,T) 2 

-A ? . 



j(q = o) 



c /-i Vskp \ mc 

c 1 1 1 v s kp \ mc 



V.c 



9v s 



node, 
antinode. 



(53) 



(54) 



We conclude that the nonlinear correction we have obtained from the q — > mode is formally two times that 
in Yip-Sauls theory. The reason for this discrepency is that we treat the response of the superconductor to A q ^ 
modes perturbatively, while the nonlinear term (v^ term in the free energy) coming from the Doppler shift is not a 
perturbative result. To see this clearly, we make a naive perturbation calculation of the response of the superconductor 
to a constant v s . The partition function is (Compare with Eq. |l4|) 



Z = exp|^Trln P (1) (iuj„) + P (2) |=cxp|^Trln P (1) {iuj n ) jexpj^TiTn 1 + pW _x (iw„)P (2) j, (55) 
where 



P«(^ n )=(- iW « k + ^ 



6k 



, P( 2 )=k 



F ■ V s 



I 
I 



(56) 



with P^ 2 ) assumed small. The free energy is 

T = -T\nZ = -T^Trln[pW(iw n )] -T^Trlnfl + pW- 1 ^)^]. (57) 

n n 

The term of present interest is expected to arise from 



rVTr 



1 



3 1 1 




PW-^JPW =-T^l f pM^|^U0. (58) 



^ 1 3 

A nonlinear term oc cannot be obtained from the perturbation theory, since the perturbation calculation implies 
v s • ki? -C Ak, while the nonlinear v% term comes from those k points near nodes on the Fermi surface satisfying 
Ak < v s • k^. Therefore the crucial step of taking (e/mc)A q ^o as v s in obtaining Eq. ( j54|) is an approximation since 
the couplings of the quasiparticles to A q ^ are treated perturbatively. Nevertheless we expect that the qualitative 
scaling behavior will be correctly reproduced with the present formalism. The renormalizations of the naive results 
to be expected from treating the magnetic field in a fully self-consistent manner are discussed in the Appendix. 

IV. PENETRATION DEPTH 

The penetration depth for a half infinite system may be defined as 

A=^y H(y)dy, (59) 

where H is the magnetic field at the surface. A specular scattering surface boundary condition on the quasiparticle 
wave functions is assumed, which replaces the surface with a current sheet of the form 

jext(y) = ~HS{y)&. (60) 

in an infinite system J2(|. 

A. Penetration depth at a constant v s 

As a first step, in this subsection, we study the penetration depth in the system characterized by an electromagnetic 
response function calculated in the presence of a constant v s . Since the true, self-consistently determined superfluid 
velocity decays with distance from the surface, this procedure introduces errors in the nonlinear terms, which we 
discuss below and estimate in some detail in the Appendix. For the moment, we consider Maxwell's Eq. 

V 2 A= (j ext +Lc) (61) 

c 

and the London Eq. (p8|), obtaining 



_ 2 f°° dq 8 f°° MC(q,v a ,r) 

Aspcc ~W 4^(q,v s ,T)/c + ^- A °~cy d<? (Ao 2 + ^2- ( 3 



Here v s is now to be understood as its value at the surface v s — v s (y = 0) = eXoH/mc as given by the solution to 
the linear, local electrodynamics problem. Obviously, A spcc includes both the nonlocal and nonlinear effects. It is 
expected to reduce exactly to the nonlocal expression of Kosztin and Leggett H if the v s dependence is neglected, 
and (qualitatively) to the nonlinear expression of Yip and Sauls [Q if the q dependence is neglected. We first go over 
these two limiting cases. 

For the linear limit v s — > 0, the qualitative behavior of the penetration depth depends on two effective nonlocal 
energy scales, E^l loc = v F u eh /X and E^ loc = v F u e i 2 /X for h,l 2 = ±1 and ug h > ugi 2 . It is shown that 

(Un) I ^( ln2 ) E I=± 1 <7T ' fOT 4olc> Clloc « T 

dX ^ ~ > i(ln2)^ (i £ + ^K-Wg,., for « T « E £l loc . (63) 

E«= ±1 «*,-««?, + 0-80 k (£) 2 E«=±i ^ iorT « E^ oc , E^ loc 
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For 8^0, E^J lloc ~ -E„onioc an d there is no intermediate parameter region ^ OI J loc <C T 1 <C £?„onioc- ^ n tms case i we 
recover Kosztin and Leggett's result [Q. However, for near any node, ^~^ loc disappears, and the linear T behavior 
of the penetration depth is preserved even at T <C vp/X , but with a reduced prefactor. We conclude that in the 
linear limit the nonlocal effects fail to cut off the linear T dependence of the penetration depth when v s is along nodal 
directions. On the other hand, it is precisely in this limit that the effects of Andreev bound states on the penetration 
depth are expected to be largest. Although the influence of these states on the field dependence has not yet been 
calculated, it appears unlikely that the naive Yip-Sauls result can apply. 
For the local limit q — * 0, Eq. ( ]6^ ) becomes 

/ \ 1/2 

A (loc) = [ | . (64) 



spcc V4^(q^0,v s ,T) 

in, for a general 9, there are two eff< 

gets 



Again, for a general 9, there are two effective nonlinear energy scales, -E^onlin = v skF u 6h an d -^nonim = v skF u eii- O nc 



, ], .,. , I 2 ( m ^) J2l=±l u 6l A ' ^nonlin' ^nonlin ^ ^ 

— ~ } l(^Ki 2 £ + 4:L « t « e„ ■ (65) 



Now we are in the position to study the penetration depth with both nonlocal and nonlinear effects. We are 
basically interested in the case of T <C v s kp. Inserting Eqs. (p3) and (Mq) into Eq. leads to 



A ^ A 



where 



with 



agL k F ul_, „ _ r 2 ^ 



771 (a;) = 4t 



fWg J -I 

—7= r^r- — -mihl) + —7= rrrr ,, 7rV2{hgi), (67) 

x 2 (3tt + 16a; 2 + 96a; 4 - 48a; 3 (l + 4a; 2 ) arctan(l/2x) 



2a; . . 

arctan(2ai) 



1 + 4x 2 



6(1 + 4a; 2 ) 



and 



47r 2 a; 4 [2 + 12x 2 -6a;(l + 4x 2 )arctan(l/2x)] 

m{x) = WT^~) ' (68) 



01 -h , h = mXoVs — -k— , (69) 



defining the competition between the nonlinear and nonlocal effects, with Hq as defined in Eq. ([5]). It is straightforward 
to get the asymptotic behavior of the penetration depth: 



(5Aspcc 7T ^ ii 1 1 1 T 2 . 

A ~ i7f "- + ce;iK 7^2 + C9i2K ^4 a 2 ' / "" L (,0) 



where the nonlinear effect dominates, with cgu = 0.0082m 3 , _ ; and cgn = 0.0059-u 3 , _ ( /itL and 

<5A^ 7r T 2 

— — ^deiiK^ 1 + dgi 2 K~~ 1 h 2 + dgizK—^, hgi < 1, (71) 

Ao 8y2 ^0 

where the nonlocal effects dominate, with dgu — O.bOu^ug^i, dgn — 1.09wg;/ l ug and dg;3 = OAlugJug.-i. In the 
local, if — > limit, the result reduces to 
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SX S pec 3 > H 

A " 2V2^W 



(72) 



which is apparently a factor of 3/2 larger than the result (jZ]). In fact, the comparison is a bit more subtle, as we have 
taken a definition of the penetration depth which differs slightly from that of Yip and Sauls. ]l2|] This comparison is 
discussed further in the Appendix. 

From Eqs. ([70|) and ( [7l| ) one now sees clearly that the linear H dependence of JAgpec at hoi ^ 1 is modified to an 
H 2 dependence when hot <C 1, implying that the nonlinear effect may be cutoff by the nonlocal effects. The crossover 
field for this to happen is defined by hgi ~ 1. For most 0, ugi and ug-i are order of unity and the crossover field is 
simply defined by ft. ~ 1, i.e., H* ~ 7TK _1 iJo/3 ~ H c i, the lower critical field. Since at H > H c i the Meissner state is 
unstable to the Abrikosov vortex state, the nonlinear Meissner effect is effectively unobservable due to the nonlocal 
effect in this geometry p2fl . In Fig. 2, we display the magnetic field dependence of the penetration depth for 6 = 0. It 
is clear that in this case the predicted linear behavior in H is only recovered well above H c \. On the other hand, for 
6 very close to a node, min(hgi) ~ 0, meaning that the crossover field can be so small that the nonlinear effect cannot 
be cut off by the nonlocal one |22| . The physical interpretation of this is that only those quasiparticle excitations near 
node regions can reduce the shielding current. When v s is along a nodal direction, those quasiparticles near this node 
acquire a large nonlinear energy shift, but on the contrary a negligible nonlocal effect since v s is perpendicular to q. 
This result suggests that the best chance for experimentalists to see the nonlinear Meissner effect is in (110) surface 
geometry. We again caution, however, that for this geometry, we do not take account of the anomalous Meissner 
currents |27j] carried by the Andreev surface bound state (2^] . In Fig. 3 we show the 9 dependence of the normalized 
penetration depth correction at T = 0. 
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FIG. 2. Magnetic field dependence of the normalized penetration depth correction SX(H,T)/Xo for 8 = 0. k= 100 is assumed 
and hence H c i/Hq = 0.01. 
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FIG. 3. 8 dependence of the normalized penetration depth correction 8X(H)/Xq for T = 0. k is the same as in Fig. 2. 



Examining the T dependence of SX in Eqs. (|7C| ) and (fn]) we find that for both h$i 3> 1 and <C 1 the leading T 
term is quadratic. In Fig. 4(a), we plot the T dependence of the penetration depth at fixed H/Hq to show explicitly 
this feature. We also display [X(H, T) — A(0, T)]/Xq in Fig. 4(b) as a function of T in comparison. Note in particular 
that the magnitude of the field dependence decreases with increasing temperature, as is to be expected for any effect 
which depends on the sharpness of the d-wave nodes. While the size of this decrease is diminished by the nonlocal 
corrections, there is never an increase in field dependence with increasing temperature, as observed in experiment (see 
discussion below). 




FIG. 4. T dependence of (a) the normalized penetration depth correction 8\{H,T)/Xq\ (b) the normalized penetration depth 
correction (X{H,T) — A(0, T))/Xo- « is the same as in Fig. 2. 



B. Penetration depth in (001) surface case 

Up to now, we have discussed in detail the penetration depth when the field is along the c-axis direction. In that 
case, all the interesting physics is in the 2D ab plane and it is implied that the possible interlayer transport kinetic 
energy is absorbed into the normal state energy spectrum. However, one may also consider the special situation in 
the quasi-2D cuprates in which the experimental configuration is characterized by a (001) surface with H || ab. In this 
case, A and v s are still in the ab plane (perpendicular to H) forming an angle 9 with the a axis. The nonlinear effect 
that has been discussed in the previous sections remains the same. However, the normal of the surface, and hence 
the penetration direction q, is perpendicular to the ab plane. If the quasiparticles are confined rigorously to the ab 
plane, they will not contribute to the nonlocal effects. But on the other hand, the magnetic field in this case is not 
screened (A — ► oo) (2jJ. In a real system, an interlayer coupling exists, leading to a nonvanishing -B„"„; oc fl^l - Here 
we confine ourselves to the situation of coherent transport along the c direction. This may be a reasonable model at 
least for the YBCO material. To avoid confusing notation, we keep k as a 2D vector in the ab plane, while we use k c 
as the momentum along the c direction. The energy spectrum of the quasiparticles in the superconducting state in 

the absence of the magnetic field is E' k kr — \/££k + ^k> where 



^ = 2^ + 2M M (73) 

with M 3> m the effective mass along the c direction. We have assumed that the pairing order parameter has no 
k c dependence. Now a parallel calculation for the penetration depth can be done. It is shown that the momentum 
q picks out k c so as to replace the first argument of the scaling function Fx in Eq. ( p3| ) by qkp c /(2MT). Thus the 
nonlocal energy scale becomes E^^ loc ~ k Fc /(MX ) ~ ^ Qc A Q /X which is no longer 9 dependent. One immediately 
sees that the main results presented for H \\ c case remain: the available nonlocal effects in (001) surface case serve 
to cut off the linear T singularity and modify the T dependence of the penetration depth to a quadratic one; at 
extremely low T, they cut off the linear H dependence when the field is below a crossover critical field H* {ab ^ which 
is found to be fp( ab ) ~ (£,Oc/£,o)H* ■ H*^ ah ^ is much smaller than H* since the c-axis coherence length is £oc — as 



13 



opposed to the in-plane coherence length of £o — 15A On the other hand, the lower critical field is also much smaller 
for this geometry, ~ 4 ^x°x @1' wnere -^Oc is the penetration depth for supercurrents along the c axis, and 

K = y/nnQ C with kq c the c-axis Ginzburg-Landau parameter. Using Ao c — 0.5 — 1 x 10 4 A, we find a large crossover 
field i/*( af) ) ~ H^ b \ making it still impossible from a practical point of view to extract a linear- If term |j3lf . 

For highly anisotropic samples like BSCCO system, the above argument might not apply. A full treatment of this 
problem awaits a generally accepted theory of the (incoherent) c-axis transport in the normal state. In addition, it 
is possible that surface-barrier effects render the field of first flux penetration much larger than we have estimated, 
leaving a large range of fields where the nonlocal effects can be neglected for this geometry. We discuss this situation 
below. 



V. DISCUSSION 

We first summarize our results. In this paper, we have investigated a clean (quasi-) 2D d x 2_ y 2-wave superconductor 
in the Meissner state based on the weak-coupling BCS theory. The existence of nodes on the Fermi surface leads to 
several important effects at low energy that have nontrivial consequences on the free energy and the penetration depth. 
In addition to the thermal excitations of quasiparticles with typical energy scale T, there are also nonlocal effects, with 
energy scale -Enonioc, due to the divergent effective size of Cooper pairs along nodal directions, and nonlinear effects 
with energy scale -E n0 niin, arising from the magnetic field-induced quasiparticle excitations. Taking the approximation 
of a slow varying supercurrent in space, we have formulated the Helmholtz free energy for the description of Bogoliubov 
quasiparticles with Doppler-shifted energy spectrum, and of the response of these quasiparticles to a weakly spatially 
varying magnetic field. The free energy is shown, after a Legendre transformation, to reduce to Schopohl-Dolgov's 
free energy at -Enoniin T, and to Volovik's singular free energy form at -Enoniin 3> T. The resulting response function 
includes both nonlocal and nonlinear effects, from which a two-parameter scaling function of the penetration depth 
5X(H; T) ~ F{E non ii n /T, E non i oc /T) is obtained. The well-known linear T dependence of 5X is obtained only at T> 
-Enoniinj -Enonioc, and will be renormalized to T 2 whenever -Enonioc or -Enoniin is larger than T. The linear H dependence 
of SX predicted by Yip and Sauls, the socalled nonlinear Meissner effect, is recovered for -Enoniin ^> T,E non i oc , but 
is typically changed to H 2 if -Enoniin is smaller than either T or E non \ oc . At extremely low T, the nonlocal effects 
cutoff the nonlinear effect at -Enonim < -E-nonioc- This happens when H < H* . Both -Enoniin and E non i oc turn out to be 
geometry dependent. When the magnetic field is along the c axis and the angle 9 that the supercurrent makes with 
the a axis is not near a nodal value, H* ~ H c i, leading to the unobservability of the nonlinear Meissner effect. When 
6 is near a nodal value, E non i oc is so small that the nonlinear Meissner effect may be recovered. However the effects 
on X(H, T) of Andreev bound states, neglected in this work, need to be better understood before this conclusion can 
be taken seriously. 

Are any of the above predictions supported by experiment? While the linear-T temperature dependence in the 
penetration depth measurements above 1 K in clean samples of cuprate superconductors has been observed for several 
years, nonlinear and nonlocal consequences of the simple d-wave theory of electromagnetic response in the Meissner 
state have proven more difficult to confirm. In zero external field, any low-T SX ~ T 2 behavior observed until now 
occurs in an experimental situation where the result can more plausibly be attributed to impurities. Recent high- 
resolution resonant coil measurements have indeed identified linear field dependences <5A ~ bH with coefficients b at 
temperatures T ~ 1 K not too different from the predictions of the local, nonlinear theory. However, the temperature 
dependence in such cases is opposite to that predicted, and the linear field dependence is not always observed. This 
suggests that earlier reports of the observation of the nonlinear Meissner effect displaying larger linear- H terms in 
lower resolution experiments must be treated with skepticism. Other manifestations of the nonlinear Meissner effect, 
e.g., the ac transverse magnetization measurements of Bhattacharya et al. pjj ] have also failed to observe the predicted 
angular harmonics, fl4| , |33|| and it is tempting to conclude that a single mechanism is responsible for the null result 
obtained by all types of experiments. 

What is the nature of this mechanism? Our proposal that the linear field dependence must generically be supressed 
at fields H ~ H& clearly does not provide a complete explanation for the failure to observe the nonlinear Meissner 
effect, as no H 2 dependence has been reported in any sample. There appear to us to be two classes of explanation for 
the existing data on 5X(H) and transverse magnetization. First, as suggested by Carrington et at, JTgj] a few vortices 
trapped in the surface layer can give rise to a "spurious" contribution to SX(H) due to the microwave excitations of 
these vortices in their pinning potentials (Campbell pinning penetration depth). ||37f Depending on the size of the 
pinning force, this can give a result, cx H, of roughly the correct magnitude. It is not clear whether the transverse 
magnetization experiments can be quantitatively explained by such a contribution. 
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The second possible explanation for the disagreement of the experiments with both the local and nonlocal d- 
wave theories was put forward by Bhattacharya et al., |2lf| who suggested the existence of a small bulk subdominant 
order parameter component, e.g., of s symmetry, which condenses in a state tt/2 out of phase with the dominant 
d component. Such a state has a true "minigap" of order the size of the s component. The authors of Ref. pi] ] 
estimated that a minigap of a few degrees K would be sufficient to eliminate the nonlinear signal in their experiments 
on YBCO. In such a scenario, any experiment which depends on the existence of quasiparticles well below the minigap 
scale must report a null result. Several other experiments have reported significant quasiparticle densities well below 
this temperature. In particular, the thermal conductivity measurements of Taillefer et al. [ j36| reveal the existence of 
a normal-fluid density in both the YBCO and BSCCO systems down to at least 50 mK. We therefore do not believe 
this proposal is tenable. 

An argument against the importance of nonlocal effects in recent experiments with H \\ ab which may provide a clue 
to the origins of the spurious signal has been given recently by Bhattacharya et al. |52| These authors point out that 
the field of first flux entry, H* x ~ 250-300 Oe, appears to be an order of magnitude larger than the Ginzburg-Landau 
lower critical field H c \ ~ 30-50 Oe, of order the crossover field between nonlocal and local field dependence for YBCO 
in this geometry. There should therefore, according to these authors, be a rather large intermediate field region before 
flux penetrates the sample where the local theory applies and the nonlinear Meissner effect should be observable. This 
argument, while intriguing, begs the question of why the observed field of first flux penetration is so different from 
the Ginzburg-Landau H c i Ipij . The obvious explanation is a geometrical surface barrier effect of the type discussed 
by Bean and Livingston, and by De Gennes. J35| ] Vortices which nucleate in the surface layer in the intermediate 
field range, while thcrmodynamically stable in the bulk, are expelled by an image potential due to the surface. The 
intermediate field "Meissner" state is therefore thermodynamically unstable, and vortices may be easily trap ped at 



the sample corners or within the skin depth. Two recent high-resolution penetration depth measurements 119 20 
failed to measure the predicted intrinsic (nonlinear) temperature dependence in this field range, but did observe a 
large linear- H field dependence most likely attributable to trapped vortices. On the other hand, the lack of hysteresis 
observed by Bhattacharya et al. when cycling above H*i argues against a significant number of pinned vortices. 
Nevertheless, the thermodynamic instability of the Meissner state suggests that the existence of a region of perfectly 
laminar surface flow with H c \ < H < H*i where the simple theories presented above apply may be too naive. The 
correct picture of the surface layer in such a situation may be one of fluctuating supercurrents whose time averaged 
effect on those quasiparticles present may be similar to that of trapped vortices. These fluctuations would make the 
field penetration layer instantaneously highly nonuniform, potentially leading to a broadening and a redistribution of 
spectral weight of the predicted transverse magnetization harmonics. fll^| We believe the only field regime where one 
can safely assume a true Meissner response is for fields below the thermodynamic H c i, i.e., in the nonlocal response 
regime. 

We propose two types of experiments which could help clarify the present impasse regarding theoretical interpre- 
tation of experimental results. The first is a direct measurement of the magnetization using a dc superconducting 
quantum interference device. This has not been yet performed as a function of applied magnetic field, to our knowl- 
edge. If the requisite sensitivity can be obtained, such an experiment would have the advantage of eliminating the 
Campbell penetration depth contributions, which may be dominating the signal in the resonant coil experiments. 
p9| , p0[ The second experiment, to maximize the possiblity of seeing nonlocal effects and eliminate the uncertainties 
in the analysis with H\\ab, would be a measurement with field H\\c on a long cylindrical sample with geometrical axis 
coinciding with the crystalline c axis. These may help to settle this long-standing controversy. 
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APPENDIX: RENORMALIZATION OF THE PENETRATION DEPTH 



In Sec. IV, we calculated the penetration depth in the presence of a constant v s . In a real system, however, the v s 
field of the screening current is varying in space. The correction to the penetration depth due to the spatially varying 
nature of v s needs to be examined in order to check the validity of the qualitative conclusions we have reached. Note 
that the rigorous result can only be obtained through the solution of the full self-consistent magnetostatic problem, 
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from which analytical information is difficult to obtain. However, we can improve our previous results by employing 
better approximations. We may imagine the system to be subdivided into many layers, in each of which v s is roughly 
constant. But v s in different layers is not the same. The response function in real space is now approximatedly defined 
as 



£(y,y') - & 



y-y ,v s 



y + y' 



(74) 



The only terms not included in this approximation involve gradients of v s . However, in the regime where E non u n 3> 
Enonioc we estimate Vv s /m <C v 2 and hence gradients of v s are negligible where the nonlinear effect is important at 
all. 

According to A = (mc/e)v s , Eq. (^TJ) can be rewritten in an operator representation 
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In the coordinate representation v s takes the form of 
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At this stage, the q and v s dependence of SIC can be approximately taken to be the same as what we got in Sec. II, 
with the point kept in mind that v s is a function of space which needs to be solved self-consistently based on Eqs. 
([78|)-(p0|). Using the definition in Eq. (|59|) and the relation H(y) = (mc/e)(dv s {y)/dy), we find 



A 



spec tt^s\ /' 

eH 



which turns out to be 
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It is easy to check that in the case of a constant r„. the penetration dept h in Eq. (N'ii mlunv exact ly n> t liat in Eq. 



(83) 



Now we insert Eqs. (33) and (f48|) into Eq. (|82|) to obtain the normalized change in the penetration depth 
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where asi = v fu$ y -i / (2^/2T) and /3#z = mv fugi / '{V2T) . Clearly the nonlinear and the nonlocal corrections to the 
penetration depth are involved in Igi(ti,R). 

For h ^> 1 it is straightforward to find that the leading correction to the penetration depth is linear in H . It is 
given by the first integral in I$i(ti, R) with v s (R) replaced by vi°\R) and the upper bound of the integral 2\hei(R)\ 
expanded to oo. After proceeding with this one gets 

SXspec Ce H -1, , ^ ,,1 /nrX 

— ^V = 2^o = J7 = 2 K K k9l>>1 ' l = ±h (85) 



Comparing with the first term in Eq. (|72|) we find that an extra prefactor 2/3 is acquired through this renormalization 
scheme. When we compare this result with Yip-Sauls result (Q), however, as we mentioned at the end of Sec. IV 
A, we note a different definition for the penetration depth is used in their paper, and an extra prefactor 2 should be 
expected. This is to say that if we use the definition of Yip-Sauls for the penetration depth, we should find it 2 times 
as large as that in Eq. (^5j). The reason that the expected 2 is actually missing here is again the consequence of the 
perturbation treatment of the coupling of the quasiparticles to the A q ^o mode, the detailed discussion of which has 
been presented in Sec. III. 

For h -C 1, the leading correction to the penetration depth is the contribution from the nonlocal effects and is 
independent of H . It can be obtained from the second integral in Iei(ti, R) by setting the lower bound 2hei(R) to be 
0. The leading independent correction to the penetration depth is found to be quadratic of H. After some algebra 
we get 

^^Vd^-i + c^VT A hm<lj i = ±h (86) 

We find that the first term, i.e., the correction from the nonlocal effects is exactly the same as the first term in Eq. 
(fn]). This means, as expected, that the nonlocal effects are not influenced by the spatial varying nature of v s . The 
second term coincides qualitatively with that in Eq. (f7l|), but acquires an insignificant extra prefactor through the 
renormalization scheme. If this expansion were continued, it is clear the result obtained would converge to the result 
of Yip and Sauls p2| in the local limit. 
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